Theory of the tunneling spectroscopy of ferromagnetic superconductors 
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We study tunneling conductance in normal metal / insulator / ferromagnetic superconductor 
junctions. The tunneling spectra show a clear difference between spin-singlet s-wave pairing, spin- 
triplet opposite spin pairing and spin-triplet equal spin pairing: These pairings exhibit, respectively, 
gap struture, double peak structure and zero bias peak in the spectra. The obtained result may 
serve as a tool for determining the pairing symmetry of ferromagnetic superconductors. 



Magnetism and superconductivity have been under in- 
tensive pursuit in the field of low temperature physics. 
Recently the interplay of them has also attracted much 
attention because nontrivial phenomena are predicted or 
found experimentally. Such phenomena are expected to 
occur in ferromanget/superconductor junctions^ 2 -! 3 - and 
also in ferromangnetic superconductors (FS) . Up to now, 
several bulk materials, e.g., UGe^, ZrZn^ and URhGc 6 , 
are identified as FS. How Cooper pairs are formed in 
FS or under the coexistence of ferromagnetism and su- 
perconductivity is an interesting problem. However the 
pairing symmetries of FS are still controversial. 

Ferromagnetic superconductors seem to be triplet su- 
perconductors because singlet pairing and ferromag- 
netism are antagonist while triplet pairing have a uni- 
form magnetic moment. However the possibility of s- 
wave pairing cannot be excluded ,L£i£ i 10 i 11 i 12 i 13 i 14 For ex- 
ample it is predicted that UGe2 can have s-wave super- 
conductivity mediated by local ferromagnetic spins i 10 ! 11 
The study of the nuclear relaxation rate cannot rule out 
the possibility of s-wave pairing in UGe2i 12113 A weak 
ferromagnetic Fermi liquid theory also suggests the pos- 
sibility of s-wave superconductivity.™ Therefore detailed 
comparison between theoretical predictions and exper- 
imantal data is required to settle this problem. Then 
the properties of thermodynamic quantities should be 
noted: For example equilibrium thermodynamic quanti- 
ties for Balian-Werthamer state of p-wave pairing, which 
is realized in B phase of 3 He, are expected to show s- 
wave property because its gap is constant^ In this way, 
equilibrium thermodynamic quantities for p-wave pairing 
could not be clearly distinguished from those of s-wave 
pairing. Therefore nonequilibrium quantities are more 
desirable to compare with experimental data. Although 
some predictions are made on the properties of junctions 
with equal spin pairing(ESP) FSf±£i±£>±&±2 the study of 
tunneling spectra for possible candidate pairings of FS is 
insufficient. 

Tunneling spectroscopy provides an important infor- 
mation on the superconducting gap and its pairing sym- 
metry. In normal metal / supercunductor junctions, 
Andreev reflection (AR)^£ is a key concept for low en- 
ergy transport. Blonder, Tinkham and Klapwijk (BTK) 
formulated the tunneling conductance where the AR is 
taken into account^!. This enables us to study the en- 
ergy gap of superconductors. The generalization of the 



BTK formula for normal metal / unconventional super- 
conductor junctions are also useful to study the prop- 
erties of unconventional superconductor o 22 ! 23 ' 24 because 
the tunneling conductance is sensitive to the pairing sym- 
metry due to the formation of midgap Andreev resonant 
states^ 3 .. 

In the present paper we study the tunneling conduc- 
tance in normal metal / insulator / ferromangnetic su- 
perconductor (N/FS) junctions. The tunneling spectra 
show a clear difference between spin-singlet s-wave pair- 
ing, spin-triplet opposite spin pairing(OSP) and spin- 
triplet equal spin pairing(ESP). This result may be use- 
ful in determining the pairing symmetry of ferromagnetic 
superconductors. 

Let us start with an effective Hamiltonian for the 
Bogoliubov-de Gennes (BdG) equation. The Hamilto- 
nian reads 
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with H (k) = + h • a , and A (k) = iAa y for singlet 
pairing or A (k) = (d (k) • a) ia y for triplet pairing. Here 
k, h and a denote electron band energy measured 
from the Fermi energy, electron momentum, applied mag- 
netic field and Pauli matrices respectively. In this paper 
we consider three types of pairings: spin-singlet s-wave 
pairing, spin-triplet OSP and spin-triplet ESP. OSP and 
ESP are characterized by the relations h x d (k) = and 
h ■ d (k) = respectively^ 

We consider a two dimensional ballistic N/FS junctions 
at zero temperature. The N/FS interface located at x = 
(along the y-axis) has an infinitely narrow insulating 
barrier described by the delta function U(x) = U8(x). 
We first consider OSP. The BdG equation reads 
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for electron-like quasiparticles, and 
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for hole- like quasiparticles, with 
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For ESP, eigenfunctions for the Hamiltonian are given by 



/ 

V 



\ 



-iO 



( \ 



/ 



u+ 














+e-» ; 







(9) 



\ 



IA 2 



\ 



1 - 



|A 2 



E±Th 



(7) 



(8) 



h = h/|h| and |A| = ±TrAA\ We assume A < h 
because otherwise the gap vanishes for the "-" state as 
can be seen in Eq. ((4|. The solution of the BdG equation 
for s-wave pairing have the same form as that of OSP and 
obtained by choosing A (k) = i Aa y . Below we consider 
unitary state for triplet superconductors and choose, as a 
model calculation, h = — hz, d (k) = A (k x + iky) /ki for 
OSP, and d (k) = A (k x + iky) /kx for ESP. Here x and 
z are unit vectors oriented to x- and z-axis respectively. 
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where is an angle with respect to the interface normal 
in the N region. Note that the magnitude of h can be 
larger than that of A for ESP because Cooper pairs are 
insensitive to the exchange field. 

We will calculate the tunnling condunctance, following 
the BTK method^ii 2 ^ Wave function ip(x) for x < (N 
region) is represented as 



ip (x < 0) 
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for an injection wave in up spin state, for s-wave pairing 
and OSP. a is AR coefficient and b is normal reflection 
(NR) coefficient. For an injection wave in down spin state 
and the junction with ESP, wave functions are given in 
a similar form. 

Similarly for x > (FS region) tJj(x) is given by the 
linear conbination of the eigenfunctions. Note that since 
the translational symmetry holds for the y-direction, the 
momenta parallel to the interface are conserved. 

The wave function follows the boundary conditions: 



Ox 



</>(+o) = v(-o), 

%(-0) = ^( + 0). 



dx 
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Applying BTK theory with AR and NR coefficients for 
electron injections with up and down spin states, we can 
calculate the angle-resolved dimensionless conductance 
for OSP represented in the form: 
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with quasiparticle energy 
E = E+ = E-, effective mass m, Fermi wavenumber 
kp and Fermi energy Ep. For s-wave pairing, the con- 
ductance is given by just replacing Y v a and Y v a with To- 
rn Ea.(fT6|). We define a^ a as the conductance in the 
normal state which is given by 

°** = 4TSJ- (19) 

The normalized conductance is represented as 
fj^ d6 cos 8 (a s+ + a S -) 
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For ESP, the conductances are given by 
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Note that (9c — \9\) have to be multiplied for a = + 

(minority spin) in Eq. (f20|) with 9c = cos -1 because 

of the mismatch of Fermi surfaces of majority and minor- 
ity spins^ Here Q(x) is the Heaviside step function. 

In the above we choose the same effective mass in N 
and FS. In most cases the effective mass in N is much 
smaller than that in FS. However it is expected that 
this effect does not change the results qualitatively for 
large Z (Z > 1)M Therefore we choose the same effective 
mass. The inclusion of the difference of effective masses 
is straightforward. 27 ' 28 Although it is known that other 
characteristics, e.g., the shape of Fermi surfaces should 
be taken into account in some phase-sensitive tests^ we 
use a cylindrical Fermi surface in this paper for simplic- 
ity because Fermi surface of FS has very complicated 
structureii&^L^ 2 . 

We study the normalized tunneling conduntace <jt as 
a function of bias voltage V . The conductances with 
Z = 10 are shown in Figs. 1(a), 1(b), and 1(c) for s- 
wave pairing, OSP and ESP respectively. For s-wave 
pairing, a gap-like sturcture appears at ft, = 0. 21 With 
the increase of h, the magnitude of the gap is reduced 
from 2 A to 2 A - 2h (Fig. Q] (a)). For OSP, a zero bias 
peak appears at h — as shown in Fig. 1(b), which stems 
from the formation of midgap Andreev resonant states^ 
We find a splitting of peak for OSP as h increases. These 
shifted structures are attributed to the h dependence of 
wave function in Eqs.([7]) and ©, and hence expected to 
emerge for all OSP (not restricted to the present choice 
of d (k)). On the other hand, the tunneling conductance 
has a zero bias peak and is almost independent of the 
exchange field for ESP as shown in Fig.l (c). This is 
because there is no energy shift in the eigenfunctions as 
shown in Eqs. (fT0|) and (fTTj) . We also find that the ab- 
sence of the shifted structure is expected for all ESP by 
calculating the eigenfunctions of the Hamiltonian with 
ESP. Therefore a clear difference between three types of 
parings can be seen. Especially when the magnitude of 
the gap A is comparable to h, the tunneling spectra are 
characterized by gap struture, double peak structure and 
zero bias peak for s-wave pairing, OSP and ESP respec- 
tively. 

A corresponding plot for Z = 1 is shown in Fig. 2. 
As shown in Fig. 2(a), the reduced dip structure ap- 
pears for s-wave pairing, the width of which is given by 
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FIG. 1: (color online) Normalized tunneling conductance with 
Z = 10 for (a)s-wave pairing, (b)OSP and (c)ESP. 



2 A — 2h. When A ~ h, the dip transforms into a single 
peak. As for OSP, a zero bias peak is formed and its 
width is reduced by the increase of h (see Fig. 2(b)). A 
zero bias peak remains with the increase of h for ESP as 
shown in Fig. 2(c). Thus there is no qualitative differ- 
ence between OSP and ESP. This is because the effect of 
midgap Andreev resonant states becomes weak for small 
Z and hence the zero bias anomaly is smeared for small 
Z. Therefore we find that the difference between for s- 
wave pairing, OSP and ESP becomes clear for large Z. 

In summary we have studied the tunneling conduc- 
tance in normal metal / insulator / ferromagnetic super- 
conductor junctions. We have found a clear difference 
in tunneling spectra between spin-singlet s-wave pairing, 
spin-triplet OSP and spin-triplet ESP. The difference is 
clear for large barrier parameter Z. This result may serve 
as a tool for determining the pairing symmetry of ferro- 
magnetic superconductors. 
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FIG. 2: (color online) Normalized tunneling conductance with 
Z = 1 for (a)s-wave pairing, (b)OSP and (c)ESP. 
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